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Magnetic fields on small scales are ubiquitous in the Universe. Although 
they can often be observed in detail, their generation mechanisms are 


not fully understood. One possibility is the so-called small-scale dynamo 
(SSD). Prevailing numerical evidence, however, appears to indicate that 
an SSD is unlikely to exist at very low magnetic Prandtl numbers (Pry) such 
as those that are present in the Sun and other cool stars. Here we have 
performed high-resolution simulations of isothermal forced turbulence 
using the lowest Pr,, values achieved so far. Contrary to earlier findings, 
the SSD not only turns out to be possible for Pr,, down to 0.0031 but also 
becomes increasingly easier to excite for Pru below about 0.05. We relate 
this behaviour to the known hydrodynamic phenomenon referred to as the 
bottleneck effect. Extrapolating our results to solar values of Pr, indicates 
that an SSD would be possible under such conditions. 


Astrophysical flows are considered to be susceptible to two types of 
dynamo instability. First, a large-scale dynamo (LSD) is excited by 
flows exhibiting helicity, or more generally, lacking mirror symme- 
try, due to rotation, shear and/or stratification. It generates coher- 
ent, dynamically relevant magnetic fields on the global scales of the 
object in question‘. The characteristics of LSDs vary depending on 
the dominating generative effects, such as differential rotation in the 
case of the Sun. Convective turbulence provides both generative and 
dissipative effects’, and their presence and astrophysical relevance 
is no longer strongly debated. 

The presence of the other type of dynamo instability, namely the 
small-scale or fluctuation dynamo (SSD), however, remains controver- 
sial in solar and stellar physics. In an SSD-active system, the magnetic 
field is generated at scales comparable to or smaller than the charac- 
teristic scales of the turbulent flow, enabled by chaotic stretching of 
field lines at high magnetic Reynolds number’. In contrast to the LSD, 
excitation of an SSD requires markedly stronger turbulence’. Further- 
more, ithas been theorized that it becomes increasingly more difficult 
to excite an SSD at very low magnetic Prandtl number Pry (refs. 4-10), 


the ratio of kinematic viscosity vand magnetic diffusivity 7. Inthe Sun, 
Pr,,can reach values as low as 10 °-10 “ (ref. 11), thus seriously repudiat- 
ing whether an SSD can at all be present. Numerical models of SSDs in 
near-surface solar convection typically operate at Pry = 1 (refs. 12-18) 
and thus circumvent the issue of low-Pr,, dynamos. 

A powerful SSD may potentially have a large impact on the dynami- 
cal processes in the Sun. It can, for example, influence the angular 
momentum transport and therefore the generation of differential rota- 
tion””°, interact with the LSD” or contribute to coronal heating via 
enhanced photospheric Poynting flux”®. Hence, itis of great importance 
to clarify whether or not an SSD can exist in the Sun. Observationally, it 
is still debated whether the small-scale magnetic field on the surface of 
the Sun has contributions from the SSD or is solely due to the tangling 
of the large-scale magnetic field by the turbulent motions” ”. How- 
ever, these studies show a slight preference of the small-scale fields 
to be cycle independent. SSDs at small Pr,, are also important for the 
interiors of planets and for liquid-metal experiments”. 

Various numerical studies have reported increasing difficulties in 
exciting the SSD when decreasing Pry (refs. 6,10,34), confirming the 
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theoretical predictions. However, current numerical models reach only 
Pry, = 0.03 using explicit physical diffusion or slightly lower (estimated) 
Pry, relying on artificial hyperdiffusion”*. To achieve even lower Pry, 
one needs to increase the grid resolution massively (see also ref. 35). 
Exciting the SSD requires a magnetic Reynolds number (Rep) typically 
larger than 100; hence, for example, Pry = 0.01implies a fluid Reynolds 
number Re = 10*, where Re = U,ms¢/V, With U,ms being the volume inte- 
grated root-mean-squared velocity, £a characteristic scale of the veloc- 
ity and Re, = PryRe. In this Article, we take this path and lower Pry 
substantially using high-resolution simulations. 


Results 

We include simulations with resolutions of 256° to 4,608? grid points 
and Re = 46 to Re = 33,000. This allows us to explore the parameter 
space from Pr,, =1to Pry, = 0.0025, which is closer to the solar value than 
has been investigated in previous studies. For eachrun, we measure the 
growth rate A of the magnetic field in its kinematic stage and determine 
whether or not an SSD is being excited. 

To afford an in-depth exploration of the effect of Pru, we omit 
large-scale effects such as stratification, rotation and shear. We avoid 
the excessive integration times, required to simulate convection, by 
driving the turbulent flow explicitly under isothermal conditions. 
Our simulation set-up consists of a fully periodic box with a random 
volume force (see Methods for details); the flow exhibits a Mach num- 
ber of around 0.08. In Fig. 1, we visualize the velocity and magnetic 
fields of one of the highest-resolution and -Reynolds-number cases. 
As might be anticipated for low-Pr,, turbulence, the flow exhibits 
much finer, fractal-like structures than the magnetic field. Note that 
all our results refer to the kinematic stage of the SSD, where the mag- 
netic field strength is far too weak to influence the flow but otherwise 
arbitrary. 


Growth rates and critical magnetic Reynolds numbers 

In Fig. 2, we visualize the growth rate À as function of Re and Rey. We 
find positive growth rates for all sets of runs with constant Pry if Rem 
is large enough. A increases always with increasing Rey as expected. 
Surprisingly, the growth rates are distinctly lower within the interval 
from Re = 2,000 to Re = 10,000 than below and above. With the Re, 
values used, this maps roughly toa Pr, interval from about 0.1 to 0.04. 

The growth rates for Pry = 0.1 match very well the ones from 
ref. 10, indicated by triangles in Fig. 2. From Fig. 2, we clearly see that 
the critical magnetic Reynolds number Rej;", defined by growth rate 
A=0, first rises as a function of Re and then falls for Re > 3 x 10° (see the 
thin black line). Looking at Rey; as a function of magnetic Prandtl 
number Pr,,, it first increases with decreasing Pr,, and then decreases 
for Pry < 0.05. Hence, an SSD is easier to excite here than for 
0.05 < Pry, < 0.1. We could even find a nearly marginal, positive growth 
rate for Pry = 0.003125. The decrease of À at low Pry, is an important 
result as the SSD was believed to be even harder*” or at least equally 
hard’* to excite when Pry was decreased further from previously inves- 
tigated values. The growth rates agree qualitatively with the earlier 
work at low Pry (refs. 6-8). : 

For a more accurate determination of Rey, we next plot the 
growth rates for fixed Pry, as a function of Rey (Fig. 3a). The data are 
consistent with A « In(Rey/Rem_ )as theoretically predicted**”. Fitting 
accordingly, we are able to determine Rey, as a function of Pry 
(Fig. 3b). This plot clearly shows that there are three distinct regions 
of dynamo excitation. When Pr,, decreases in the range 1 > Pry, > 0.1 it 
becomes much harder to excite the SSD. In the range 0.1 > Pry, 2 0.04, 
excitation is most difficult with little variation of Ref". For Pry < 0.04, 
it again becomes easier as Pr, reduces. In refs. 7,8, the authors already 
found an indication of Reĝ" to level-off with decreasing Pry, however, 
only when using artificial hyperdiffusion. Similarly, with our error bars, 
a constant Rey; cannot be excluded for 0.01 < Pry, < 0.1. However, at 


Pry, = 0.005, the error bar allows to conclude that Regit is here lower 


than at Pry = 0.05. This again confirms our result that Regi" is decreasing 


with Pry for very low Pry. : 

For Pry < 0.05, the decrease of Rey, with Pry can be well repre- 
sented by the power law Regi" x Pry”. Extrapolating this to the Sun 
and solar-like stars would lead to Rey; ~ 40at Pry = 10%, which means 
that we could expect an SSD to be present. For increasing Re, by 
decreasing v, it would be reasonable to assert that the statistical proper- 
ties of the flow and hence Rey become independent of Pry. However, 
episodes of non-monotonic behaviour of Rey; when approaching this 
limit cannot be ruled out. l 

The well-determined Rey; dependency on Pry together with 
its error bars and the power-law fit have been added to Fig. 2, and 
agree very well with the thin black line for A = 0 interpolated from the 


growth rates. 


Regions of dynamo excitation 

Next we seek answers to the obvious question arising: why is the SSD 
harder to excite in a certain intermediate range of Pry and easier at lower 
and higher values? For this, we investigate the kinetic and magnetic 
energy spectra of arepresentative subset of the runs (Supplementary 
Table 2). We showin Fig. 4 the spectra of two exemplary cases: run FOOS, 
with Pr,, = 0.05, probes the Pr,, interval of impeded dynamo action, 
while run HOOOS, with Pry = 0.005, is clearly outside it (see Supple- 
mentary Figs. 1 and 2 for spectra of other cases). 

In all cases, the kinetic energy as a function of wavenumber k 
clearly follows a Kolmogorov cascade with Fj, < K°” in the inertial 
range. When compensating with k*”, we find the well-known bottleneck 
effect”: alocal increase in spectral energy, deviating from the power 
law, as found both in fluid experiments*°’ and numerical studies****. 
It has been postulated to be detrimental to SSD growth”. For the 
magnetic spectrum, however, yet clearly visible for only Pry < 0.005, 
we find a power law following Emag * K°. A3/2 slope at low wavenumbers 
as predicted by ref. 45 is seen only in the runs with Pr,, close to one, 
while for the intermediate and low-Pr,, runs, the positive-slope part of 
the spectrum shrinks to cover only the lowest k values, and the steep 
negative slopes at high k values become prominent. A steep negative 
slope in the magnetic power spectra was also seen by ref. 7 for Pry 
slightly below unity. However, the authors propose a tentative power 
of—1given that the -3 slopeis not yet clearly visible for their Pry, values. 

Analysing our simulations, we adopt the following strategy. For 
each spectrum, we determine the wavenumber of the bottleneck, k,, 
as the location of its maximum in the (smoothed) compensated spec- 
trum, along with its starting point ky, < k, at the location with 75% of the 
maximum (Fig. 4, middle). We additionally calculate a characteristic 
magnetic wavenumber, defined as ky = JiEmag(K)KAK/SEmag(k)dk, which 
is often connected with the energy-carrying scale. Furthermore, we 
calculate the viscous dissipation wavenumber k, = (ex/v3) “following 
Kolmogorov theory, where ex is the viscous dissipation rate 2vS? with 
the traceless rate-of-strain tensor of the flow, S. From the relations 
between these four wavenumbers (listed in Supplementary Table 2), 
we draw insights about the observed behaviour of Rey, with respect 
to Pry. 

We plot k,/k, and k,,/k, as functions of Pr,, in Fig. 5. As is expected, 
k,/k,, or the ratio of the viscous scale to the scale of the bottleneck, 
does not depend on Pry, as the bottleneck is a purely hydrodynamic 
phenomenon. The start of the bottleneck k,, should likewise not 
depend on Pry, but the low Re values for Pry =1 to Pry = 0.1 lead to 
apparent thinner bottlenecks, hence an unsystematic weak depend- 
ency. The red shaded area between k, and ky, is the low-wavenumber 
part of the bottleneck where the slope of the spectrum is larger (less 
negative) than -5/3 (see Supplementary Table 2 for values of the modi- 
fied slope a, and Supplementary Section 1 for a discussion). We note 
that a, = -1.3 ... -1.5 and can thus deviate markedly from —5/3. Overplot- 
ting the k,,/k, curve reveals that it intersects with the red shaded area 
exactly where the dynamo is hardest to excite (region II). This lets us 
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Fig. 1| Visualization of flow and SSD solution. Flow speed (left) and magnetic field strength (right) froma high-resolution SSD-active run with Re = 18,200 and 


Pry, = 0.01 0n the surface of the simulation box. 


rail i COTENT 
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Re 


Fig. 2 | SSD growth rate as function of the fluid and magnetic Reynolds 
numbers (Re and Re,,). The diamonds represent the results of this work and the 
triangles represent the results of ref. 10. The colour coding indicates the value of 
the normalized growth rate Àt with T= 1/U,msKp a rough estimate for the turnover 
time. The dotted lines indicate constant magnetic Prandtl number Pry. The white 
circles indicate zero growth rate for certain Pry, obtained from fitting for the 
critical magnetic Reynolds number, as shown in Fig. 3; fitting errors are signified 
by yellow-black bars (Supplementary Section 5). The background colours, 
including the thin black line (zero growth), are assigned via linear interpolation of 
the simulation data. The green dashed line shows the power-law fit of the critical 
Rey for Pry < 0.08, with power 0.125 (Fig. 3b). 


conclude that the shallower slope of the low-wavenumber part of the 
bottleneck may indeed be responsible for enhancing Rey, in the inter- 
val 0.04 < Pry, < 0.1. Using this plot, we can now clearly explain the three 
regions of dynamo excitation. For 0.1 < Pry < 1 the low-wavenumber 
part of the bottleneck and the characteristic magnetic scale are com- 
pletely decoupled. This makes the SSD easy to excite (region I). For 
0.04 < Pry, < 0.1, (grey, region II), the dynamo is hardest to excite 
because of the shallower slope of the kinetic spectra. In region III, where 
Pry < 0.04 the low-wavenumber part of the bottleneck and the charac- 
teristic magnetic scale are again completely decoupled making the 
dynamo easier to excite. 

Further, we find that the dependence of k,,/k, on Pry also differs 
between the regions. In region I, ky/k, depends on Pry via ky/ky « Prov 


and in region lI and III via ky/k, « Prei ' This becomes particularly 


interesting when comparing the characteristic magnetic wavenumber 
ky with the ohmic dissipation wavenumber which is defined as 
k, = k Prii. Inregion I, we find a notable difference of kyand k, in value 
and scaling. However, in region III, the scaling of ky comes very close 
to the 3/4 scaling of k,. This behaviour can be even better seen in the 
inset of Fig. 5, where the ratio k,,/k, is 0.3 for Pry = 1and tends towards 
unity for decreasing Pry, but is likely to saturate below 0.75. 


Discussion 

In conclusion, we find that the SSD is progressively easier to excite for 
magnetic Prandtl numbers below 0.04, in contrast to earlier findings, 
and thus is very likely to exist in the Sun and other cool stars. Provided 
saturation at sufficiently high levels, the SSD has been proposed to 
strongly influence the dynamics of solar-like stars: previous numerical 
studies, albeit at Pr, = 1, indicate that this influence concerns, for exam- 
ple, the angular momentum transport” and the LSD”. Our kinematic 
study, however, only shows that a positive growth rate is possible at 
very low Pry, but not whether an SSD is able to generate dynamically 
important field strengths. As the Rey of the Sun and solar-like stars is 
several orders of magnitude higher than the extrapolated Rey, value 
of 40, we yet expect dynamically important SSDs as indicated by Pry, =1 
simulations. However, numerical simulations with Pry down to 0.01 
show a decrease of the saturation strength with decreasing Pry, (ref. 46). 

The results of our study are well in agreement with previous 
numerical studies considering partly overlapping Pry ranges® *”°. 
Those studies found some discrepancies with the Kazantsev theory* 
for low Pry, for example, the narrowing down of the positive Kazantsev 
spectrum at low and intermediate wavenumbers, and the emergence 
of a negative slope instead at large wavenumbers’. We could extend 
this regime to even lower Pr,,and therefore study these discrepancies 
further. For Pry < 0.005, we find that the magnetic spectrum shows a 
power-law scaling K?°, which is substantially steeper than the tentative 
Kone proposed in ref. 7 for 0.03 < Pry < 0.07 (but only for eighth-order 
hyperdiffusivity). This finding of sucha steep power law in the magnetic 
spectrum challenges the current theoretical predictions and might 
indicate that the SSD operating at low Pry is fundamentally different 
from that at Pry = 1. 

Second, we find that the growth rates near the onset follow an 
In(Rey,) dependence as predicted by refs. 36,37, and nota Rei? oneas 
would result from intertial-range-driven SSDs’”. We do not observe a 
tendency of the growth rate to become independent of Re,, at the 
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Fig. 3 | Growth rate and critical Reynolds number. a, Normalized growth rate At 
as function of magnetic Reynolds number Rey for simulation sets with fixed 
magnetic Prandtl number Pry, indicated by different colours. Logarithmic 
functions At « In(Rey/Rexr) according to refs. 36,37 were fitted separately to 
the individual sets, as indicated by the coloured lines (see the dashed-dotted line 


for the mean slope). b, Critical magnetic Reynolds number Reg" as function of 


150 L 


crit 
M 


100 L 


Re 


sob 


0.10 
Pr, 


1.00 0.01 


M 


Pry obtained from the fits in a. The error bars show the fitting error 
(Supplementary Section 5). The diamond indicates a run with growth rate À = 0; 
hence, its Rey represents ~ Regi" for the used Pry = 0.003125. The red dashed line 
is a power-law fit Regt" x Pro, valid for Pry, < 0.08. The grey shaded area 
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indicates the Pry, interval where the dynamo is hardest to excite (Rey 2 150). 
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Fig. 4 | Energy spectra. Kinetic (top) and magnetic (bottom) energy spectra for 
two exemplary runs with Re = 7,958 and Pry = 0.05 (left), and Re = 32,930 and 

Pru = 0.005 (right). In the middle row, the kinetic spectra are compensated by kK°®?. 
Vertical lines indicate the forcing wavenumber k; (green solid), the wavenumber 
of the bottleneck’s peak k, (red solid) and its starting point k,, (red dotted), the 
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viscous dissipation wavenumber k, (orange), the ohmic dissipation wavenumber 
ky = kPt (dark blue) and the characteristic magnetic wavenumber ky (light 
blue). All spectra are averaged over the kinematic phase whereupon each 
individual magnetic spectrum was normalized by its maximum, thus taking out 
the exponential growth. 


highest Pr,, either, which could be an indication of an outer-scale driven 
SSD, as postulated by ref. 7. Furthermore, we find that the pre-factor 
of y x In(Rey/Rey) is nearly constant with its mean around 0.022, in 
agreement with 0.023 of ref. 10. A constant value means that the loga- 
rithmic scaling is independent of Pr,, and seems to be of general 
validity. 

Third, we find that the measured characteristic magnetic wave- 
number ky is always smaller than the estimated k,, and furthermore, 
ky does not always follow the theory-predicted scaling of k, œ Pret 
with Pry. For region I, where Pr,, is close to 1, this discrepancy is up toa 
factor of three and the deviation from the expected Pr,, scaling is most 
pronounced here. These discrepancies have been associated with the 
numerical set-ups injecting energy at a forcing scale far larger than the 


dissipation scale, that is k; « k, (ref. 1). Furthermore, our runs in region 
l also have relatively low Re and therefore numerical effects are not 
dismissible. In region III (low Pry), ku/k, is approaching the constant 
offset factor 0.75. Hence, the scaling of ky/k, with Pr, gets close to the 
expected one. This result again indicates that the SSD at low Pr, is dif- 
ferent from that at Pr, ~ 1. 

Anincrease of Rey; with decreasing Pry followed by an asymptotic 
levelling-off for Pry « 1was expected in the light of theory and previous 
numerical studies. Instead, we found non-monotonic behaviour as 
function of Pry; we could relate it to the hydrodynamical phenomenon 
of the bottleneck. If the characteristic magnetic wavenumber lies in 
the positive-gradient part of the compensated spectrum, where the 
spectral slope is markedly reduced from -5/3 to about -1.4, the dynamo 
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Fig. 5 | Relation of the characteristic magnetic wavenumber k,, to the 
bottleneck. We show its peak k, and its starting point k,, in red, the characteristic 
magnetic wavenumber kyin light blue and the ohmic dissipation wavenumber k, 
in dark blue. The red shaded area between k, and k,, corresponds to the 
low-wavenumber part of the bottleneck where the turbulent flow is rougher than 
for a -5/3 power law. The Roman numbers indicate the three distinct regions of 
dynamo excitation. The region of the weakest growth (II) is over-plotted in grey. 
The characteristic magnetic wavenumber ky, can be fitted by two power laws 
(black dotted lines): ky/k, « Prov for Pry 2 0.05 and ky/k, « Pro for 

Pry < 0.05. All wavenumbers are normalized by the viscous one k,. We find that 
the dynamois hardest to excite if kų lies within the low-wavenumber side of the 
bottleneck. Leaving this region towards lower or higher wavenumbers makes the 
dynamo easier to excite. The inset shows k,,/k, as a function of Pry. 


is hardest to excite (0.1 > Pry 2 0.04). For higher or lower Pry, the 
dynamo becomes increasingly easier to excite. The local change in 
slope due to the bottleneck has often been related to an increase of the 
‘roughness’ of the flow'°, which is expected to harden dynamo excita- 
tion based on theoretical predictions*” from kinematic Kazantsev 
theory*. Inline with theory, the roughness-increasing part of the bot- 
tleneck appears decisive in our results, however, only when k,, is used 
as a criterion. The usage of k, would in contrast suggest that the peak 
of the bottleneck is decisive”. Such interpretation appears incorrect, 
as the rough estimate of k, employed here does not represent the 
magnetic spectrum adequately and the peak of the bottleneck does 
not coincide with the maximum of ‘roughness’. 


Methods 

Numerical set-up 

For our simulations, we use a cubic Cartesian box with edge length L 
and solve the isothermal magnetohydrodynamic equations without 
gravity, similar to refs. 5,47. 


Du 


Tia -c2VInp +J x B/p + V - (2pvS)/p + f, (1) 
ðA 
—= B+ 7V2A, 2 
aE uxB+n (2) 
Dp 
Dim -V - (pu), (3) 


where u is the flow speed, c, is the sound speed, p is the mass density, 
B = V x Ais the magnetic field with A being the vector potential and V 
is the gradient vector. J = V x B/u, is the current density with magnetic 
vacuum permeability 4o, while vand 7 are constant kinematic viscosity 


and magnetic diffusivity, respectively. The rate-of-strain tensor S;= 
(uij + u;;)/2 - 6,;V - u/3 is traceless, where ô; denotes the Kronecker 
delta, and the Einstein notation convection applying to their indices 
i and j. The forcing function f provides random white-in-time 
non-helical transversal plane waves, which are added in each time step 
tothe momentum equation (see ref. 5 for details). The wavenumbers 
of the forcing lie in a narrow band around k;= 2k, with k, = 211/L. Its 
amplitude is chosen such that the Mach number Ma = u,,,./C, is always 
around 0.082, where u,ms = ,/(u2), is the volume and time-averaged 
root-mean-square value. The Ma values of all runs are listed in Sup- 
plementary Table 1. To normalize the growth rate À, we use an esti- 
mated turnover time T = 1/(Ujmsk; )= 6/(k,C,). The boundary conditions 
are periodic for all quantities and we initialize the magnetic field with 
weak Gaussian noise. 

Diffusion is controlled by the prescribed parameters v and 7. 
Accordingly, we define the fluid and magnetic Reynolds numbers 
with the forcing wavenumber k; as 


Re = Upms/Vkp, Rep = Urms/Nke. (4) 


We performed numerical free decay experiments (Supplementary 
Section 7), from which we confirm that the numerical diffusivities are 
negligible. 

The spectral kinetic and magnetic energy densities are defined via 


| Exin(k) dk = u? ps ()y/2, (5) 
k 


I Emag(k) dk = Brms/2Ho. (6) 
k 


where Bims = \/ (B°), is the volume-averaged root-mean-square value 
and (p),is the volume-averaged density. 

Our numerical set-up employs a markedly simplified model of tur- 
bulence compared with the actual one in the Sun. There, turbulence is 
driven by stratified rotating convection being of course neither isother- 
mal nor isotropic. However, these simplifications were so far necessary 
when performing a parameter study at such high resolutions as we do. 
Nevertheless, we can connect our study to solar parameters in terms 
of Pr,,and Ma. Their chosen values best represent the weakly stratified 
layers within the bulk of the solar convection zone, where Pry < land 
Ma < 1. The anisotropy in the flow on small scales is much weaker there 
than near the surface and therefore close to our simplified set-up. 


Data availability 

Data for reproducing Figs. 2, 3 and 5 are included in the article and 
its Supplementary Information files. The raw data (time series, spec- 
tra, slices and snapshots) are provided through IDA/Fairdata service 
hosted at CSC, Finland, under https://doi.org/10.23729/206af669- 
07fd-4a30-9968-b4ded5003014. From the raw data, Figs. 1 and 4 can 
be reproduced. 


Code availability 

We use the Pencil Code** to perform all simulations, with parallelized 
fast-Fourier-transforms to calculate the spectra on the fly”. Pencil Code 
is freely available at https://github.com/pencil-code/. 
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